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Abstract: Searches for lepton flavour violating r — > 111 and r — > I7 decays at the 
B-factories are presented. Upper limits on the branching ratios are obtained ~ 
O(10 -7 ) at 90% confidence level. 



1. Introduction 

Lack of any pre-existing experimental evidence led the Standard Model 
(SM) to be built upon the assumption that the lepton flavour number is 
conserved within each family. Recent data on neutrino oscillations [f ] have 
started to provide hints of lepton flavour (LF) violation in the neutral lep- 
ton sector. Various extensions of the SM to include non-zero neutrino mass 
predict manifestations of charged LF violation in terms of branching ratio's 
of t — > 111 (1 7 ) decays at the level of fCT 14 (40 > [2]. However, many beyond 
the SM processes, eg. MSSM + heavy Majorana neutrino with seesaw mech- 
anism [3], MSSM with soft SUSY breaking terms [4], extra gauge bosons 
(technicolor) [5] predict LF violation at the level of O(10 -7 — 1 -10 ) [6]. 

Present generation e + e _ B-factories also serve as r factories owing to 
the large production cross-section of t + t~ pairs at center-of-mass (CM) 
energy around the T(4S) resonance. Recent results from the r — > 111 
search [7] in all six possible charge conserving decay channels (e - e + e~, 
/i + e~e~, [j,~e + e~ , e + ^ _ /i _ , e"/^ + ^~ and with 91.6 fb _1 data col- 

lected with the BABAR experiment is described in detail here. Results from 
the r — > ^7 search [8] with 86.3 fb _1 data collected with Belle experiment 
are also quoted. 



*To appear in the proceedings of Lake Louise Winter Institute 2004 on Fundamental 
Interactions (LLWI 2004), Lake Louise, Alberta, Canada, 15-21 Feb 2004. 
^On behalf of the BaBar Collaboration. 
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2. r — ► 111 search 

Event Reconstruction: LF violating r decays are generated according 
to flat phase space distribution, while the other tau of t+t~ pair produced 
in c+c~ annihilation decays according to the known branching fractions as 
modeled by KK2F generator including radiative effects [9]. 

Events are required to have 4 well reconstructed charged tracks, with 
each of the 3 tracks (from the signal r decay) separated by > 90° in the 
CM frame from the remaining track (from the other r decay). Pairs of 
oppositely charged tracks identified as photon conversions in the detector 
material with an e + e _ invariant mass < 30 MeV are vetoed. About 50% of 
the reconstructed MC signal events pass this 1-3 topology requirement. a 

Each event is then classified into 6 different decay modes of the signal r 
based upon its charge, the ratio of calorimeter energy to track momentum, 
the ionization loss in the tracking system and the shape of the shower in 
the calorimeter (for electron identification), the number of hits in the muon 
chamber and energy deposits in the calorimeter (for muon identification). 
The electron and muon identification have an efficiency per lepton of 91% 
and 63% respectively for the signal MC. The mis- identification rate for a 
hadron selected as an electron or muon in SM 3- prong tau decays is 2.2% 
or 4.8% respectively. 

Background Suppression: The dominant backgrounds for e~e+e~, 
e~^i + /i~ channels are Bhabhas, for [i~e + e~ , channels are di-muons, 

and for ,u + e~e~, e + n~[i~ channels are hadronic events (qq) with the extra 
tracks coming from un-identified single or multiple conversions. 

The QED backgrounds (Bhabha, di-muons) are suppressed by requiring 
the CM momentum of the 1-prong track < 4.8 GeV. For e~e + e~ and 
e~/i + /i~ channels, the 1-prong track is required not to be an electron, while 
for /j,~e + e~ and channels it must not be a muon. For all these 4 

channels, the CM angle 6*13 between the 1-prong momentum and the vector 
sum of the 3-prong momenta must satisfy cos #13 > —0.9999, and the net 
transverse momentum of the four tracks > 100 MeV. 

To reduce qq and SM t + t~ backgrounds, events in the four channels 
specified above are required to have no unassociated calorimeter clusters 
(photons) in the 3-prong hemisphere with energy > 100 MeV in the labo- 
ratory frame, while events in all six channels are required to have no track 
in the 3-prong hemisphere that is also consistent with being a kaon. 



a This estimate includes the 85% branching fraction for the 1-prong r decay. 
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Signal Yield estimation: The distinguishing feature of r — > 111 search 
is that the signal events cluster around (0,0) in (M rec - m T , E* cc - E£ oam ) 
= (AM, AE) plane, where M rec is the reconstructed invariant mass of the 
three tracks, m r = 1.777 GeV is the tau mass [10], E* oc is the total CM 
energy of the three tracks and E£ cam is the beam CM energy. Detector 
effects as well as initial and final state photon radiation broaden signal 
distributions toward lower values of AE and AM respectively as shown 
in Figure 1. Signal yield is estimated inside the smaller rectangular boxes 
shown as inset in the same figure, chosen so as to give the smallest expected 
upper limits on the branching fractions in the background-only hypothesis. 

The background events have distinctive distributions in the (AM, AE) 
plane: qq events tend to populate the plane uniformly, while QED back- 
grounds are restricted to a narrow band at positive values of AE (owing 
to the presence of extra charged particles in the event), and t+t~ back- 
grounds are restricted to negative values of both AE and AM (as the signal 
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Figure 1 : Observed data as dots and 
the boundaries of the signal region 
for each decay mode. The dark and 
light shading indicates contours con- 
taining 50% and 90% of the selected 
MC signal events, respectively. 
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Table 1: Efficiency estimates, num- 
ber of expected background events, 
and their relative systematic uncer- 
tainties, number of observed events, 
and branching fraction upper limits 
for each decay mode. 
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topology reconstruction does not account for the missing neutrino(s)). 

The expected background rates for each channel are determined by fit- 
ting a set of probability density functions (PDFs) to the observed data in 
the (AM, AE) plane in a grand sideband (GS) region, defined as the outer 
boundary of the regions shown in Figure 1 (excluding the signal region). 
The analytical forms of the PDFs used are described in detail in Refer- 
ence 7. For both qq and t + t~ backgrounds, the shapes of these PDFs are 
determined by fits to qq and t+t - background MC samples for each decay 
mode. For the QED backgrounds, the shapes of the PDFs are obtained by 
fitting data control samples having a 1-3 topology (without any restriction 
on cos #13), which are enhanced in Bhabha or events by requiring 

that the particle in the 1-prong hemisphere is identified as an electron or 
muon. 

With the shapes of the three background PDFs determined, an un- 
binncd maximum likelihood fit to the data in the GS region is performed 
to obtain the expected number of background events (-/Vb g d) in the signal 
region, shown along with selection efficiencies (e) and relative systematic 
uncertainties on e (Sys e ) and on iVbgd (Sys bgd ) in Table 1. The leading 
contributions to Sys £ come from statistical precision of data control sam- 
ples used in particle identification and modeling of tracking efficiency, while 
estimate of Sys bgd is limited by finite data statistics in the GS region used 
to estimate iVbgd- 

No significant excess is found in any decay mode, and in Table 1 
are shown upper limits on the branching fractions calculated as B^j° L — 
N^l/(2eLa TT ), where is the 90% CL upper limit for the number of 
signal events when iV bs events are observed with A^, g d background events 
expected for L = 91.6 fb _1 , a TT = 0.89 nb b including all uncertain- 
ties [11] [12]. 

3. Summary 

LF violating r decays have not been observed as yet. Upper limits at 90% 
confidence level on the branching ratios are obtained: 

B?° l (t ^ 111) ~ 1 - 3 x 10" 7 [7], Bfjl(r - 1 7 ) = 3.1 x HT 7 [8]. 
These limits represent an order of magnitude improvement over the previous 
experimental bounds. 



b The luminosity is measured using the observed (J, + [i rate, the cross-sections are esti- 
mated using KK2F [9] and the uncertainty on the product L cr TT is 2.3%. 
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Experimental limits are starting to probe predictions from several theo- 
retical models with B ~ O(10~ 7 — 10~ 10 ). Improvements over the last few 
decades on experimental limits on are shown in Figure 2 along with 
some theoretical predictions from different LF violating models. 
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Figure 2: Evolution of experimental bounds (Bfj L ) and some predictions. 
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